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We have measured the photon statistics of pump and probe beams after interaction with Rb atoms 
in a situation of Electromagnetically Induced Transparency. Both fields present super-poissonian 
statistics and their intensities become correlated, in good qualitative agreement with theoretical pre- 
dictions in which both fields are treated quantum-mechanically. The intensity correlations measured 
are a first step towards the observation of entanglement between the fields. 



PACS numbers: 42.50.-p, 42.50.Ar, 42.50.Ct, 42.50.Gy 



Electromagnetically Induced Transparency (EIT) jjj is 
an interference effect that can be observed when three- 
level atoms interact simultaneously with two lasers. The 
recent interest in these systems has been driven by ob- 
servations of very slow light pulse propagation in EIT 
media |^-|| and even of light storage ||U ■ Even though 
it is normally viewed as a quantum interference effect, 
EIT has a very simple classical counterpart The 
question of whether there exist intrinsic quantum ef- 
fects (with no classical analog) in EIT naturally arises 
in this context. It is even more important in connection 
with suggested applications in the field of quantum in- 
formation P,[h1. Quantum field properties in EIT have 
been theoretically investigated. Phase-noise squeezing 
was predicted by Fleischhauer and co-workers jllj . Agar- 
wal |l^] found matched photon-statistics for two classi- 
cal and two quantum fields interacting with three-level 
atoms, in a situation of Coherent Population Trapping. 
Jain [ p"3| extended this work, and predicted excess-noise 
correlations in EIT. In all the above p|-|i~3|], at least one 
of the fields was treated classically. As we will see below, 
it is our belief that such an approximation is not valid in 
the EIT situation. 

In this Letter we present the first experimental in- 
vestigation (to our knowledge) of field fluctuations and 
correlations in EIT. We have performed photon statis- 
tics measurements of both pump and probe fields, as a 
function of the probe field detuning. As predicted by 
our theoretical treatment, both fields, treated quantum- 
mechanically, have super-poissonian statistics. Further- 
more, they are coupled by their interaction with the 
atoms and thus their intensities become correlated. 

We begin by describing our theoretical model for three- 
level atoms in the A-configuration (ground states |1) and 
|2), and excited state |0}) interacting with two quantum 
modes of the electromagnetic field. In order to simu- 
late the interaction of the atoms with two propagating 
fields, we use the input-output formalism |l4||l5| ], and 
consider the interaction between the two fields and the 
atoms in a ring cavity E^ ], with external input fields. 



For the intracavity field operators we used plane wave 
and quasi- monochromatic approximations. The interac- 
tion hamiltonian is obtained with the usual dipole and 
rotating-wave approximations 



h gi §f(t)Ai(t) + h 92 s+(t)A 2 {t) + h.c. 



(i) 



where A± (A 2 ) is the annihilation operator for intracavity 
field 1 - pump (intracavity field 2 - probe), g\ (g 2 ) is 
the atom - field 1 (field 2) coupling strength, and 
(S 2 ) the atomic polarization on the transition |1) «-> |0) 
(|2)H0». 

From the Heisenberg equations of motion, we derive 
quantum Langevin equations fll7|| f or the system opera- 
tors. In a matrix form, we find JL8| 



dX(t) 
dt 



= -A(t) X(t)+T>f-{t) 



(2) 



The vector operator X(t) has as its components the op- 
erators giving the atomic inversions and polarizations 
(and their hermitian conjugates) corresponding to the 
two transitions, the coherences between the two ground 
states and the annihilation and creation operators of both 
fields. The vector of the Langevin forces acting on the 
system is given by J~(t), and D is the diffusion matrix. 
The matrix A(i), in the steady state, will represent the 
drift matrix. Eq. (|^) is a compact form of writing twelve 
coupled differential equations for the operators. 

We are interested in studying the fluctuations of the 
field operators. This is done by linearizing the operators 
in Eq. (||) around their stationary values. A new set of 
twelve coupled differential equations is obtained. It is in- 
teresting to examine the equations for the fluctuations of 
one field and for the corresponding atomic polarization: 



dSA 1 
dt ' 

dSSi 



iA cl 5 Ax 



i — db 1 
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SA U 



(3) 



i8li ) 5S 1 + igiwiSAx + igiaiSWi 

(4) 



dt V 4 

-ig 2 s* 12 SA 2 - ig 2 a 2 5S^ 2 + F S i 
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Here we define 7 as the cavity linewidth, A c i cavity de- 
tuning for field 1, r cavity length divided by the speed of 
light, Ann the annihilation operator for the input field 1, 
Ti spontaneous emission rate from |0) — > |1), Sli detun- 
ing between field 1 and the corresponding atomic tran- 
sition, wi steady-state inversion between states |0) and 
|1), ai (012) steady-state amplitude of field 1 (field 2), 

inversion (operator) between states |0) and |1), s 12 
steady-state coherence between ground states |1) and |2), 

coherence operator, F$i Langevin fluctuation force. 
The notation JS-f means fluctuations of the correspond- 
ing operator. 

It is clear from Eqs. (|J) and (Q) that the fluctuations 
in field 1 are not only determined by input field fluctua- 
tions but also by atomic fluctuations and by fluctuations 
in field 2. Moreover, in Eq. ([!]) we notice that noise cor- 
relations between the fields will arise, conditioned to the 
existence of a coherence between the ground states, s* 2 , 
only important on EIT resonance. 

The terms igiWiSAi, igiaiSWi, «<72Si2<^2, ig2Ct2<5S^2 
in Eq. (|J) arc all of the same order. This means that 
we can not neglect the fluctuations of either field with 
respect to its average (steady-state value) . It is not valid 
to treat one field as classical and the other as quantized. 
This conclusion also applies to the work by Agarwal |l2] ] 
and Jain Jl3| who consider four independent fields (two 
classical and two quantized) interacting with coherent 
population trapped (CPT) atoms. The fields used to 
produce CPT also have inherent quantum fluctuations, 
which lead to atomic fluctuations and, therefore, should 
not be neglected. 

From the twelve coupled differential equations for 
the fluctuations, we calculate the spectral matrix 
S(fl) jl9|,^0| for the noise, which yields the fluctuations 
of the output cavity fields, after Fourier transform of the 
input-output relations. 

The use of a cavity to simulate an EIT experiment with 
cw propagating fields (no cavity) is adequate as long as 
we restrict ourselves to fluctuations around the steady 
state and do not investigate transient effects. We must 
also consider the cavity large enough so that the atoms 
experience neither significant changes in their sponta- 
neous emission rates nor collective effects. 

The noise spectra calculated for probe and pump am- 
plitude quadratures and their correlations are plotted in 
Fig. [l], as a function of the probe frequency. These spec- 
tra were calculated for an analysis frequency f2 = r/67r 
(consistent with the experiment), where T is the to- 
tal spontaneous emission rate from the excited state. 
We took equal coupling constants for both transitions, 
pump strictly on resonance, and a ratio of intensities 
[|ai/a2| 2 = 9 in Fig. |l] (a)] such that there is a very 
narrow and deep EIT resonance. For the correlations in 
Fig. [I] (b), both fields have equal intensities. The num- 
ber of atoms was of the order of 10 4 , which is sufficiently 



large to justify the linearization method used. 

We find super-poissonian statistics for both fields, 
and intensity correlations peaked at the EIT resonance. 
Physically we understand this behavior as follows. The 
coupling of the fields with the atoms introduces a spread- 
ing in the photon distributions associated to each field. 
The atoms act as "beamsplitters", redistributing pho- 
tons between both modes, in such a way that the mean 
numbers do not vary appreciably but the variances are 
modified. Taking absorption into account, we expect this 
effect to be reduced, since absorption introduces random- 
ness in the photon redistribution process. We therefore 
expect to observe a maximum effect for minimum absorp- 
tion, which corresponds to the EIT resonance. The pre- 
dicted correlations also depend on the presence of absorp- 
tion. In the EIT situation, the transmission of one field 
depends on the (fluctuating) intensity of the other field 
(cross Kerr effect), leading to intensity correlations. The 
experimental results presented below are in very good 
qualitative agreement with these predictions. 
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FIG. 1. Theoretical predictions, (a) Super-poissonian pho- 
ton statistics for both fields, as a function of probe detuning. 
The left vertical axis is for the pump (gray) and right axis for 
the probe (black), (b) Intensity correlations as a function of 
probe detuning (in units of V). 
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FIG. 2. Sketch of the experimental setup. DL1 and DL2: 
extended-cavity diode lasers; FR: Faraday rotator; GP: glass 
plate; PBS: polarizing beamsplitter cube; S.A.: spectrum an- 
alyzer. 

Our experimental setup is sketched in Fig. |2[ Measure- 
ments were performed in a 5 cm long room-temperature 
Rb vapor cell. Two independent extended-cavity diode 
lasers (ECDLs) provide the pump and probe fields. A 
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small portion of each beam is extracted and sent to aux- 
iliary saturated absorption cells, used as frequency refer- 
ences. The pump beam is tuned to the 85 Rb |5Si/ 2 , F = 
3) — > I5P3/2, F' — 3) transition, while the probe beam is 
scanned across the Doppler-broadened \5Si/ 2 ,F = 2) — > 
|5P 3 / 2 , F' — 1, 2, or 3} line. The two beams have orthog- 
onal polarizations and are combined by means of a po- 
larizing beamsplitter cube (PBS). At the output of the 
cell we can separate the beams again with another PBS. 
We then have the option of detecting either beam with 
a balanced detection setup |21|. Half- wave plates also 
enable us to recombine the two fields and measure the 
photon statistics of the sum and difference intensities of 
the two beams. From these we can extract the intensity 
correlations. 

The usual EIT signals arc observed by sending only 
the probe beam into the detection region and measuring 
the DC (average) intensity. The intensity transmitted 
through the cell has a narrow peak as a function of the 
probe frequency. The width is much narrower than the 
Dopplcr width, which is a signature of the interference 
effect. Next, we measure the photon statistics (inten- 
sity fluctuations) for the probe beam, yielding the sig- 
nal presented in Fig. || (a). For the range of frequencies 
spanned, the photon statistics is always super-poissonian, 
but it presents a sharp peak (almost 20 dB) correspond- 
ing to the EIT resonance. In this measurement, the ini- 
tial photon statistics of the probe beam is poissonian (for 
the frequency of analysis chosen, the pump initially has 
super-poissonian statistics, but the theory predicts this 
effect even when both fields have poissonian statistics). 
The pump photon statistics presents similar behavior, 
as seen in Fig. || (b). In Fig. || (c) and (d) we present 
the corresponding Fano factors for both probe and pump 
fields, respectively. The Fano factor is given by the ra- 
tio of the intensity fluctuations to the average intensity 
(shot noise level), and we plot it on a linear scale. In 
these measurements, the probe and pump intensities were 
14.5 mW/cm 2 and 63.8 mW/cm 2 , respectively. The ob- 
served behavior agrees very well with the theoretical pre- 
dictions of Fig. 0. 

The noise spectra were recorded using a spectrum an- 
alyzer (HP 8560 E) in the zero span mode (RBW = 
300 kHz, VBW = 3 kHz), with a center (analysis) fre- 
quency of 2 MHz. This frequency is in a window such 
that the probe laser has poissonian photon statistics and 
the electronic noise is safely below shot noise (more than 
5 dB). The pump laser's noise is ~ 8 dB above shot 
noise (and shows little variation with the analysis fre- 
quency). The analyzer is triggered by the same signal 
used to scan the probe frequency. The photodetectors 
are EG&G FND-100 with a nominal efficiency of ~ 70%. 
We checked that none of the measurements is influenced 
by saturation in the detection. This is done by introduc- 
ing neutral density filters before detection and observing 
that the photon statistics becomes poissonian, with a lin- 



ear dependence on absorption. 

For typical EIT signals, in the presence of a strong 
pump field, the shot noise level also peaks, as expected, 
following the increasing DC intensity on EIT resonance. 
However, if we lower the pump power, this effect tends 
to disappear, since incoherent effects then dominate over 
the coherent pumping on the transition \5Si/ 2 , F = 3) — > 
|5-f3/2, -F' = 3). The probe DC signal, and consequently 
the shot noise level, no longer show any evidence of a 
coherent effect in the atomic medium for pump intensi- 
ties lower than 0.45 mW/cm 2 , for a probe intensity of 
13.3 mW/cm 2 (it is the ratio of pump to probe intensi- 
ties that is relevant). On the other hand, by looking at 
the photon statistics, the coherent effect is still clearly 
identifiable. 
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FIG. 3. (a) and (b) Noise spectra of the probe and pump 
beams, respectively, as a function of the probe frequency 
and corresponding shot (in gray) and electronic noise (lower 
traces), (c) and (d) Fano factors deduced from (a) and (b). 

We have also measured the noise properties of the sum 
and difference beam intensities. This is done by send- 
ing each beam to one of the detectors. By subtracting 
the intensity-difference noise from the intensity-sum, we 
observe a clear correlation on the EIT resonance. The 
corresponding shot-noise level is measured by mixing the 
two beams so that each detector receives half of each 
beam. The results are presented in Fig. ||. We observe 
~ 7 dB splitting between the sum and difference fluc- 
tuations. For photon statistics, this corresponds to a 
correlation term 2C = 4{(mri2) — (ni)(ri2)}. 
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FIG. 4. (a) Measurements of the noise in the sum (black) 
and difference (gray) intensities of pump and probe beams. 
The difference observed on EIT resonance is a measurement of 
the intensity correlation (b), on a linear scale, created between 
both fields. Pump and probe intensities: 7.8 mW/cm 2 and 
7.5 mW/cm 2 , respectively. 

In summary, we have measured the photon statistics of 
probe and pump beams after interaction with an atomic 
medium, in a situation of Electromagnetically Induced 
Transparency. We observe super-poissonian statistics, 
peaked on EIT resonance. We have also observed inten- 
sity correlations between the two initially independent 
fields, as a result of their interaction with the atoms. 
The experimental result features agree very well with our 
theoretical predictions. The balanced detection scheme 
leads to enhanced sensitivity for detecting coherent ef- 
fects in the atomic medium. The importance of the cor- 
relations predicted and observed is two-fold. On the one 
hand, this is a very strong indication that, when investi- 
gating quantum properties of one of the fields in a coher- 
ent situation (such as EIT), one can not treat either field 
as classical, as done by several authors. The coherent 
exchange of photons between the fields, which is at the 
very origin of the EIT effect, introduces correlations be- 
tween their fluctuations. Another important issue is the 
nature of the correlations between both fields. Quan- 
tum correlations (entanglement) between intense fields 
can be used in the context of quantum information. This 
would be the dual of the experiment by Julsgaard and 
co-workers |22| , in which correlations between atoms in 
two different vapor cells were created using light. In the 
present experiment, we have made the first step in this 
direction by demonstrating the existence of correlations. 
However, in order to distinguish between classical and 
quantum correlations, we will need an improved exper- 



imental setup. A new calculation shows that entangle- 
ment between pump and probe fields can be produced in 
the EIT situation, and its observation will require mea- 
surements of quadrature fluctuations of both fields |^3| . 
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